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DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


OFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIE 
N ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH TH 
ION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER 
Pp 
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Ooo 
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$ 
; 
S_THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
' tt al NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 


THE INFORMATION IN THIS SOFTWARE IS SUBJEC 
AND SHOULD NOT BE CONSTRUED AS A COMMI 
CORPORATION. 


DIGITAL ASSUMES NO RESPONSIBI 
SOFTWARE ON EQUIPMENT WHICH I 


fofololojle=lo) 
Sooooo 
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TMENT BY DIGITAL EQUIPMENT 
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S NOT SUPPLIED BY DIGITAL. 
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; FACILITY: MATH LIBRARY 

; ABSTRACT: 

This module contains the routine MTHSGMOD: 

It returns the remainder of the division of argl/arg2 using 


the following equation: 
argl - Cint(argl/arg2))*arg2 


8S SS SS SSS SSS SSSSSSSSSSSsesss 
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; AUTHOR: Jeffrey C. Wiener, CREATION DATE: 21-DEC-1982 
; MODIFIED BY: 
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.SBTTL DECLARATIONS 

: INCLUDE FELES: 

NONE 


EXTERNAL SYMBOLS: 


DSABL GBL : Force all external symbols to be dectared 
“EXTRN MTHSSSIGNAL 

SEXTRN MTHSK_FLOUNDMAT 

SEXTRN MTHSK"INVARGMAT 

! LIBRARY MACROS CALLS: 

‘ SSFDEF : Define SF$ (stack frame) symbols 
: EQUATED SYMBOLS: 


EXP 52 
HIGA_MASK 


: OWN STORAGE: 

NONE 
; PSECT DECLARATIONS: 

~-PSECT _MTHSCODE PIC, SHR, LONG, EXE, NOWRT 
; CONSTANTS: 


TWO_EXP_52: 
~-LONG “%X00004350, “x0 3 2e#52 


*x000 


: 52%2°4 
“XFFF 


"oO 
oo 
“Ww 
nS 
"oO 


[=] 

[=] 

So 

[=] 
SoOoOSSCOSCOSOSOSOSSOOOSOOCOOSOOSOOOSOSoOO 


COONAN WN SO OONOAULS WIN (OC OWONAUE WO DONO 
. 


©96909 09 69 09 C9 09 C9 09 SI INI NINN IO A AAAAAAAO 


SOoOooooooooooooooooSao 


SOoOOooooooooooooooo 
CGOoooooooooooooooo 


oo 


00000000 00004350 


SOSOSSSOSCSOSOOSCSOSOSOSOOOSOSOSOOSOOOSOSOSOOSOSSOSOSSOOSOSOSOOOOOOOOSOSOOoOO 
0.00 09.60.09 0d Cd Cd Cd Cd 09 G9 0D CD Cd Cod 09 CD Od G9 OD GD Gd 09 G9 OD OD Od CD Cd GD CD CD G9 CD CD CD GD CD OD CD C9 C9 C9 CO. 0D 0D 09 09 09 09 C0 CD. CD OD0D0D 


SOSSSSOSSOOSOSCOSSOOSOOOOOSOOOOCSOSOOSOSOOOOSOOSOOOOOOOOO 


SoOOoSoOoooooooooooooooooooooooooooooooooooo 


MTHSGMOD - G REAL*8 remainder 6-SEP-1 


i 


NOUS WN SO OD NOAUE WIN @ O OD NAME WIN 3 O OO IA UE WIN) SO CONAUE WN OONOUI 


Be Ge Se Ge Se Ge Ge Be Se Se Se Se Sete 


FERRERS SE PUWWWWWWIWIWIWIRIRINININININININNY BQO OOOO OOOOO OOOO 


a a ek a a a a td a to = to = = = = = Ss 2 2) sd —) 2s —) —» 3d 2 a — 2 ss ss 


BSP 1984 41:58:86 EMTMRTL. SReoMTHGROD MAR; 1 


4 
‘ »SBTTL MTHSGMOD - G REAL*8 remainder 
FUNCTIONAL DESCRIPTION: 


Return the remainder of orgi/ergs in G  satcte, point format 


Remainder = argi - (int(argl/arg2))*arg 


The algorithm used to evaluate the GMOD function is as follows: 


X = the first argument. 
Y = the second argument. 
step 1. m = the exponent of Y. 
n = the exponent of X. 
c=n-™m 
If ¢ < 0, end with result = X. 
step 2. 1 = the fractional part of X. 
J = the fractional part of Y. 
If lor J, Lele d 
Go to — . 
step 3. L = 2*(p-1)*1, where p = 53 for G_floating numbers. 
step 4. T = L/J 
} 2 Pree Spee wertrs T is int(L/J) or int(L/J)+1 
s = ® 
If 1<0,1leIl¢J T was int(L/J)+1 
step 5. c = c¢ = (p=1) 
If ¢ > 0 go to step 3. 
step 6. If c_= =-(p-1) go to step 9. 
step 7. L = 2*(p-l4c) * 
step 8. l=L-J#T 
step 9. Result = 2*m * | 


CALLING SEQUENCE: 


Remainder .wg.v = MTHSGMOD (dividend.rg.r, divisor.rg.r) 


INPUT PARAMETERS: 


The two input parameters are G_floating-point values. 
passed by reference. 

DIVIDEND = 4 : Dividend = 
DIVISOR 8 : Divisor = 


IMPLICIT INPUTS: 
NONE 
FUNCTION VALUE: 


Remainder of the division argi/arg2 is returned as a 
G_floating point value. 


IMPLICIT OUTPUTS: 
NONE 


; COMPLETION CODES: 


X in the algorithm. 
Y in the algorithm. 


Sete te Se 


on the top of 


7E 52 6€ 43FD SUBG3 (SP), R2, (SP) (SP) = J2@=J- J) 


52 50 S1FD CMPG }=—s RO, R If I<J 
19 BLSS STEP go to STEP_S 
50 52 42FD SUBG2 R2, R else I = =) 
6D 14 BGTR STEP_ go to STEP_5 if 1>0. or 


else the algorithm ends 
TSTW sory teemmtars 


be the sign of the result is 
MNEGG RO, RO 


the same as the sign of 
the first argument, A. 


MTHSGMOD 16-SEP-1984 01:29:21 AX/VMS Macro v04-00 P 
3-001 MTHSGMOD = G REAL*8 remainder patie Pi 933338) MTHRTL.SRCIMTHGMOD.MAR; 1 si dh 
148 ; 
| 143 ; NONE 
13 ; SIDE EFFECTS: 
138 ; MTHS_INVARGMAT = Invalid argument to math Library if the divisor is zero. 
154 ; MTHS _FLOUNDMAT = Ang hg underflow in math Library is signaled if 
132 : the FU bit is set in the callers PSL. 
hi 
O1FC o ; 133 ~ENTRY MTHSGMOD, “M<R2, R3, R4, RS, R6, R7, RB 
5208 BF 50te bor ee <TR 2 te Sy sO. bivioten b 
: Y20. 
50 04 BC SOFD 0011 188 MOVG  @DIVIDEND(AP), RO ion <t... 
56 52 FFFF8OOF 8F CB 0016 165 BICL #*°XFFFFSOOF, R2, R6 3; R6=m is the biased t of Y 
5850 FFFFGOOF BF CB OOIE 166 BICL3 #*XFFFFBOOF, RO, RB iE ibex Sa Gee Blened eee ae 
Ss SCE 0 168 SUBL2 R6, R8 3 R4 = ¢ = n-m unbiased 
01 rf 60 $ 198 BGEQ STEP_2 3; Result is xX if x<Y, ie, if c<O 
04 4 ? 109 R 3; RO/R1 = X 
56 DD 09 ; 76 STEP_2: PUSHL R6 : push m onto the stack 
52 FFFO 8F AA 09 4 174 BICW2 #*XFFFO, R2 ; R2/R3 = J = unbiased ifract(Y): 
52 4000 8F AC O33 10? XORW #*x4000, R2 3; J = properly biased ifract(Y): 
50 FFFO 8F AA ists 177 BICW2 #*xXFFFO, RO ; RO/R1 = I = unbiased ‘fract (x): 
50 4000 8F AC 3878 178 XORW #*x4000, RO ; I = properly biased ‘fract(Xx): 
Bs 
04 182 ; In STEP_4 and STEP_8 the calculation of I = L - J#int(L/J) must be 
094 184 : computed as precis ly as possible. To do this we will need to write J as 
be ! 5 ; where J1 = the high 26 bits of J and J2 = J - J1, the low 27 bits of J. 
$b 1 ; : HIGH _MASK is used to extract the 5 bits of J from longword2 that belong 
004 188 ; to JT. 
One 189 ; 
5 
43 52 7D 7 136 MOVQ R2, -(SP) (SP) = J 
04 AE FFFFO7FF 8F CA 323 32 BICL #HIGH_MASK, 4(SP) (SP) = J1 replaces the value 
4D 195 
a 
19 
5 138 
; i 
1 


a a Be 
50 50 52FD 006 
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3-001 MTHSGMOD - G REAL*8 remainder pa oe ts 7 9:33:81 MTHRTL.S Re SRTH THGMOD .MAR; 1 = (3) 
04 es 5 DONE: RET 
79 5 ERROR: ASHQ #15 ; Y=0. Reserved operand. 
43 00'8F 9A MOVZBL #MTH an ARGHAT -(SP) : error code 
00000000 ' GF 1 re cays #1, G*ATHSSSIGNAL 3 signal the error 
50 00000340 8F CO STEP_3: ADDL2 #EXP_52, RO ; RO/R1 = L = 2e#(p-1)*] 


TEP_4: 
ee = 2*(52) is added and then subtracted from 
T = int(L/J) to ensure that T = chopped(L/J) or chopped(L/J)+1 


Se Ge Ge Ge Ge Ge Ge 


56 50 52 47FD DIVG3 R2, RO 3 R6/R7 = T = L/J 
56 OFF 77 CF 40rD ADDG2 TWO_EXP_ gS : RO/R7 = T = T+2*8(p-1) 
56 FF71 CF 42FD SUBG Two> EXP” de 3 Te2**(p-1) = yer chopped 


The calculation of I = L = J*int(L/J) must be grammed as precisely 
as possible. To ,99 this we will need to write T 


T= 
where 21 = the high 26 bits of T and 72 = T = 21, the low 27 bits of T. 
Now, using J = J1 + J2, 


OOOO oCooooooooooooooooooooooooo 


SOGOCOOCOOSOOCOSOOSOSOOOOOOOOOOOOOoOoOo 
>>> > > > >> > > se 


WOOOOOCOSOOOCCOCOCSOOCCOO®S LT - “1 “CDOS 
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009 L=-J * imt(L/J) = Lb = (1 + U2) © (21 + 22) 
009 =L = (21 * J1) = (21 * JQ) 
009 - (22 * 31) _ (22 * 92) 
09 =L-(271"))2>(2"5) 
09 
09 - 
009 
54 DO 009 MOVL. R6, R4 ; 
55 57 FFFFO7FF 8F CB 009 BICL3 #HIGH MASK, R7, RS > RG/RS = 71 
5654 42rd 98 SUBG ,R : R6/R7 = 22 
7E 08 AE 34 45FD MULGS R4, B(SP), -(SP) : Compute 21*J1 
5 E 42FD SUBG2 (S$P)+, RO > RO/RI = L - gist 
54 3 6E 44FD MULG (SP) R4 + R4/RS = Z18J 
50.54 42rD SUBG R : RO/RI = L = 21*J 
54 56 08 AE 45FD MULG BSP) R6, R4 t R4/RS = 22eJ 
50 54 42rD SUBG2 R4, RO : RO/RI=L - gisJ - 22*J1 
36 gE 44FD MULG ah R6 t R6/R7 = 22*J 
0 42FD SUBG Rg : RO/RI = L = 2 
06 14 BGTR STEP 
8 13 BEQL ng tune : If RO/R1=0 the algorithm ends 
50 52 40FD ADDG2 R2, RO 3; Add J back in because you had 
$ T=chopped(L/J)+1 
58 00000340 8F c STEP_5: SUBL2 #EXP_52, RB sc = co(p-1) = c-52 
A BGEG STEP 
58 00000340 8F C ADDL2 #EXP_52, RB zc = (p-l)4c = S52ec 
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The next two Lines of code are STEP_6 and STEP_7. 


BEQL STEP_9 $ 
ADDL2 RB, RO : L = 182*(c+t) 


2*(p-1) = 2*°(55) is added and then subtrac 1 i from 
T = int(L/J) to ensure that T = chopped(L/J) or chopped(L/J)+1 


DIVG3 R2 3 R6/R7 = T = 

ADDG2 Two Eb 8S : R6/R? =T= Fegesipa1) 

SUBG TWO_EXP_52, Re 3 ee “1) =t/Jc opped 
3 /J chopped + 

STEP_8: 


The calculation of I = L = J#int(L/J) must be ranean as precisely 
as possible. To 499 this we will need to write T 


Z 
where Z1 = “the high 26 bits of T and 72 = T = 21, the low 27 bits of T. 
Now, using J = J1 + J2, 


L=-J * int(L/J) = Lb = (1 + J2) © (21 + 22) 
=L = (21 J1) = (21 * 4g) 
- (72 * J1) | (23 * J2) 
=L-(21* J) = (22" 9) 
MOVL. 6, R4 
BICL3 #HIGH_MASK, R7, R5 : R4/RS = 71 
SUBG r + R6/R? = 722 
MULG3 R4, B(SP), -(SP) : Compute 21*J1 
SUBG2 (SP)+, RO + RO/R1 = L = Z21*J1 
MULG2 (SP), R4 t R4/RS = Z18J 
SUBG2 4 + RO/RI =L - disJ 
LG3 © B($P), R6, R4 : RG/RS = 22eJ 
SUBG2 R4, RO : RO/RI=L - gies - 22*J1 
MULG2 (SP), R6 > R6/R7 = Z2eJ 
SUBG2 6, RO : RO/RI = L - Z*J 
BGTR STEP 9 
BEQL sy : If RO/R1=0 the algorithm ends 
ADDG2 R2, RO 3; Add J back in because you had 
3 T=chopped(L/J)+1 
_9: SUBW #14000. 16(SP) ; Remove the bias from m and 
ADDW pA od: RO 5 form RO/R1 = 2°m*L 
BLSS NDER FLOW : Branch if oaineaiiaes 


TEST_SIGN: 


TSTW @DIVIDEND (AP) ; the sign of the result is 


"SSEP=1984 91:58:82 EATHNTL. seeamrncnopsmar;1 "28° 8) 


J 16 
TH 16-SEP-1984 01:29: AX/VMS M v04-00 Pp 
a ste MTHSGMOD = G REAL*8 remainder gr SEF= 1382 9] :S3:8) FATA NS acre ee obomar:1 29 (fy 
2% 1 BGEQ RETURN 3; the same as the sign of 
50 8000 8? A 133 BISw2 #°x8000, RO ; the first argument. XN. 
04 1? RETURN: RET 
14 UNDERFLOW: 
50 7C 014 4 CLRQ RO ; Set default result to 0.0 
0 04 AD 06 ET 9145 $25 BC ASFSV_FU, SFSW_SAVE_PSuCFPD, NOfUL * rs a Bs A bar? 
é ranc caller has n n 
90000000'8F pp ol4a 359 PUSHL  #MTHSK_FLOUNDMAT : Report MTHS FLOUNDMAT. 
00000000'GF 01 FB 129 CALLS #1, G*ATHSSSIGNAL i Signal the Condition 
04 01 NO_FU: RET 
1 4 0 
1 1 END 


MTHSGMOD 
Symbol table 
DIVIDEND 000004 
DIVISOR = $ 
DONE R 0s 
, 2 * 

HIGA_ MASK = FFFFO/FF 
THES 1GNAL eeeneree = =6X OO 

HS$GMOD 00000008 “ 0 
ATHSe _FLOUNDMAT ‘teeeeeee 0 
MTHSK— _INVARGMAT teeeeeee . 8 
NO_FU- 157 R 
RETURN : ie i 0 
SFSV_FU = 0 
SFSW_SAVE_PSW = i444 4H 
STEP_ 00000: Ns R 0 
STEP_ 00000078 R 0 
STEP_ 0Q000000CB R 0 
STEP_ 000001 f R 0 
TEST SIGN 00000138 R 0 
TWwO_ExP_52 44 60 R 02 
UNDERF LOG 00000143 R 02 
PSECT name Allocation 
» 6.x 00000000 ( 0.) 

$ 00000000 ( 0.) 
MTHSCODE 00000158 <(« 344.) 
Phase Page faults CPU Time 
Initialization 9 0:00:00.11 
Command processing 103 0:00: 3 
Pass 116 _ 701.34 
Symbol table sort 0 :00: 8.06 
Pass 2 71 :00:00. 
Symbol table output 3 8 
Psect synopsis output :00. 
Cross-reference output : 8-8 
Assembler run totals 326 :02.81 
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Attributes 
NOPIC 
NOPIC 


PIC 


hold 48 non-local and 
object records in Pass 2. 
macros. 


"SE 


USR 
USR 
USR 


pages) of sivtael’ eoaaré were used to buffer the intermediate 
0 pages of symbol table rt dt 
roducin 


EP-19 
P=1984 


CON 
CON 
CON 


ode. 
local symbols. 
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-ABS 
ABS 
REL 


1: 
1: 
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LCL 
LCL 
LCL 


SHR 
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NOSHR NOEXE NORD NOWRT NOVE 
NOSHR EXE RD WRT NOVE 
RD NOWRT NOVE 


ree 
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! Macro Library statistics ! 


eeeeeoeeoooooooccooecescoos} 
Macro library name Macros defined 

_$255$DUAZ8: [SYSLIBISTARLET .MLB; 2 sgeummerees e: 

88 GETS were required to define 4 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DISABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHGMOD/OBJ=OBJ$:MTHGMOD MSRC$:MTHGMOD/UPDATE=(ENHS:MTHGMOD) | 


NT CORPORATION 
D PROPRIETARY 
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